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The last decade has seen a strong upsurge in the study of quantum mechanical transport 
processes, especially based on extensions of their classical counterparts JH|2l|3l|4l|5l. On one 
hand this is due to the growing attention being given to quantum information and to potential 
quantum computers [6|; on the other hand, the interest is also fuelled by recent experimental 
breakthroughs in the coherent energy transfer over atomic and molecular systems (especially 
at extremely low temperatures) [HI, processes which, say in the framework of Frenkel 
excitons, are very closely related to classical random walk models 0. 
Previous work ifTOl fTTT fT2l IT3l has highlighted the close relations between the classical 
continuous-time random walks (CTRWs) lfT4l 1731 [Toll and the continuous-time quantum 
walks (CTQWs). One difference between CTRWs and CTQWs lies the fact that quantum- 
mechanically time is connected to the imaginary constant. Hence, the basic physical aspects 
of the classical and the quantum picture differ, although with regards to the spatial coordinates 
one finds in both cases (discrete) versions of the Laplace operator. Hence, many of the 
mathematical tools used to handle Laplacian forms in classical physics may well be used 
to bear fruit in the quantum mechanical environment. Furthermore, for simple underlying 
geometries the quantum problem is directly related to well-known models in polymer and in 
solid state physics. For instance, walks on (one-dimensional) chains are readily treated by 
a Bloch ansatz ifTTl [P71 fT8ll . when periodic boundary conditions (PBCs) are implemented. 
On the other hand, open boundary conditions (OBCs) appear naturally in polymer physics in 
relation to the Rouse-model lfl9l . 

In this work we focus on CTQWs over two-dimensional, finite networks. These are 
topologically quite simple systems which, however, can display complex quantum mechanical 
features. By keeping the systems as simple as possible, we highlight the complex behaviour 
of the quantum mechanical transport compared to the classical one. Now, assuming either 
PBCs or OBCs (separately for each direction) leads to structures topologically equivalent to 
toroidal, cylindrical, and finite rectangular networks. For such networks the eigenvalues of the 
system can be obtained in simple, closed form. The average transition probabilities between 
distinct sites display quite unexpected, odd behaviours in the long-time limit for some special 
network sizes. These findings extend the results previously obtained for square networks lfT2l 
to rectangular ones. As an additional feature, we analyse the average probability to be still 
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or again at the initial site; in this case one has an (analytical) lower bound, that in many 
ways is quite close to the numerically established behaviour. The advantage of the lower 
bound expression is that it depends only on the eigenvalues, but not on the eigenvectors of the 
system. 

The paper is structured as follows: In the next section we recall the general properties of 
CTQWs on networks. In Sec.|3]we analyse the role of the PBCs and OBCs on the transport. 
Section |4| is devoted to the average probability to be still or again at the initial site, both in 
the classical and in the quantum case. In Sec. |3 we analyse the long time behaviour of the 
transport between pairs of sites, first for square networks, and then for general, rectangular 
networks. We summarise the obtained results in Sec. HO 

2. Continuous time quantum walks on graphs 

In this work we focus on coherent quantum mechanical transport over systems whose 
topology can be modelled by two-dimensional graphs. Specifically, we will consider networks 
which consist of M x N = J\f nodes jointed by identical bonds; such networks are 
topologically equivalent to finite, rectangular lattices whose length in the ^-direction is M, 
and in the y-direction N. Then we denote the position of node j by (j x ,jy), where j x and 
j y are integer labels in the two directions. We will consider different situations: in each 
direction we either impose PBCs or OBCs. In two dimensions this leads then to three distinct 
topological objects, namely to a rectangle, to a cylinder and to a torus, see Fig.[T] Furthermore, 
we will assume that the excitation is initially localised at site j. Classically, CTRWs are 
described by the master equation ll20ll2T1 : 



where Ph,j(t) is the conditional probability to find the walker at time t at node k when starting 
at time at node j. We assume an unbiased CTRW such that the transmission rates 7 of all 
bonds are equal. Then the transfer matrix of the walk, T = (Tkj), is related to the connectivity 
matrix by T = —7 A, where A = (A k j) is a discrete form of the Laplacian operator. One 
has namely Aij = — 1 if nodes i and j are connected by a bond and A+j = otherwise; 
furthermore An = — Aij or equivalently An = fi where fi is the functionality of site 



The quantum mechanical extension of a CTRW, the CTQW, is now defined by identifying 




(1) 



1. 
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the Hamiltonian of the system with the (classical) transfer operator, H = — T HI ITUll. 
A convenient way to implement this idea is to start from localised states, so that state 
\j) = \j x ) (g> \ j y ) = \j x ,jy) is localised at node j. The states \ j) span the whole accessible 
Hilbert space and form an orthonormal basis set, i.e., (k\j) = 5kj and J2j \ = 1» with 
1 being the identity operator. If at the initial time to = only state \j) is populated, then 
the transition amplitude (Xk,j{t) from state \j) to state \k) in time t obeys the Schrodinger 
Equation (SE) 

^"fcjW = Y^ H kiUi,j(t), (2) 

whose solution can be formally expressed as: 

a h j(t) = (k\exp(-iHt)\j). (3) 

The corresponding transition probability during time t is nk,j(t) — \ a k,j(t)\ 2 - Nevertheless, 
despite the formal similarity between Eqs. CO) and 0, we have classically ^ k Pk,j{t) = 1 but 
quantum-mechanically J2k \ a k,j(t)\ 2 = 1- 




Figure 1. Two-dimensional networks with (a) finite rectangular, (b) cylindrical and (c) 
toroidal topology. These drawings do not exhaust the broad panoply of geometrically possible 
realizations of such networks and are meant solely to display clearly the connectivities. 



We hasten to note that there is no unique way of defining a CTQW. For example, for regular 
networks, where all nodes have the same functionality, different choices of the Hamiltonian 
can give rise to the same quantum dynamics [4]. Nevertheless, in what follows we will stick 
to directly identifying the Hamiltonian with the transfer operator, since some of the networks 
we will consider below are not regular. For instance, the functionality of sites of a finite 
rectangular network ranges from 2 for a corner site to 4 for an internal site. Sites on the 
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boundary of a cylinder have functionality 3. A torus is regular in the sense of Ref. 0, since 
all sites have functionality 4. 

It is now reasonable to work in an orthonormal basis \q n ), n G [1,A/] which diagonalises 
A, and hence also H = — T = 7 A. For this we need to know all the eigenvalues and the 
eigenvectors of A. As is well-known, the matrix A is non-negative definite; moreover, for 
a connected underlying structure it has only one vanishing eigenvalue, the other eigenvalues 
being strictly positive ll22l . Thus A = QAQ -1 , where Q is the matrix built up by the \q n ), 
Q : _1 is its inverse, and A is the diagonal matrix having as diagonal elements the eigenvalues 
A n of A. 

Let us now turn our attention to the probability to be still or again at the starting point of the 
motion Pk,k(t)- We can write the formal solution of Eq. (HJ in Dirac notation as 

Pk^t) = (fc|exp(Ti)|j>, (4) 

where classically the vector \ j) has a 1 at the position corresponding to the pair (j x ,j y ) and 
zeros otherwise. Then, Pk,k{t) reads 

PhAt) = ( k \ exp(T*)|fc) = ^(fc|<z n > exp(- 7 *A n )(qr n |fe). (5) 

n 

As also shown in Il23ll24ll25l . the average of the classical probability Ph,k{t) over a ^ s i tes °f 
the graph is 

P(t) = Jf^PkM 

k 

= j^^2^2(qn\k)(k\q n )exp(-jt\ n ) = ^ exp(-7A n t). (6) 

n k n 

The classical p{t) decays monotonically from p(t) = 1 to a final, asymptotic plateau, 
lim t ^oop(t) = 1/jV. This behaviour is typical for a diffusive process leading to energy 
equipartition. From Eq. © it is evident thatp(t) depends only on the eigenvalues and not on 
the eigenvectors of T. 

For CTQWs, the average probability ft (t) to be still or again at the initial state at time t reads, 



k k 



2 



^^exp[-i7(A n - X m )t\^2\(k\q n ) | 2 | (k\q ri 



(7) 
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where we used that a kjk (t) = ^2 n (k\q n ) (q n \k) exp(— i^\ n t). Eq. © depends explicitly on 
the eigenvectors, which renders its determination cumbersome. As shown in [ 13 1, 7f (t) admits 
a lower bound, as follows using the Cauchy-Schwarz inequality: 



i i M 



E 



> 



M (t). 



(8) 



Furthermore, in analogy to Eq. © we have 



k n 



^2exp{-i^X n t)(q n \q n ) = ^ exp(-«7A n t). 



(9) 



Thus fx(t) equals 



1 

A72 



Y^kAt) 

k 

^^exp[-i7(A n - \ m )t]. 



^exp(-i7A„t) 



(10) 



The last expression shows that also fi(t) depends only on the eigenvalues. As we will discuss 
in Sec. El fi(t) provides important information about the network, despite of being only a 
lower bound . 

3. Role of different boundary conditions 

We now turn to two-dimensional networks built from J\f = M x N nodes, which, by 
implementing different boundary conditions (OBCs or PBCs) stay rectangular or turn into 
cylinders or tori, see Figd These symmetries are based on the form of the matrix A and are 
only topological in nature; in fact, it is not compulsory for these structures to obey any kind 
of translation symmetry, see for instance the situation in polymer physics [26| or in quantum 
chemistry ll27l . where the A-matrices appear in Hiickel molecular orbital calculations 0. 
For the sake of simplicity and without any loss of generality we set now the transmission rate 
7 = 1. 

To fix the ideas let us first consider the one-dimensional case, namely a finite chain consisting 
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of M nodes, arranged in the x-direction. Under PBCs the corresponding Hamiltonian, H£, 
reads 

ni\j x ) = 2\j x )-\j x -i)-\j x + i), (ii) 

where this relation is to be understood modulo M, in particular \M) = |0) and \M + 1) = 
|1). The corresponding eigenfunctions are \^e x ) = (1/VM) Y^f x =i ex P( — i9 x j x )\j x ) and 
correspond to the eigenvalues \ dx = 2 — 2cos^ x , where 9 X = 2imi/M, with m = 
0,1, ...,M- 1. 

Under OBCs the corresponding hamiltonian H° reads: 

H°|l> = |1) - |2> , H° X \M) = \M) - \M - 1) (12) 

and 

n° x \j x ) = ni\j x ) (13) 

for \j x ) with j a . G {2, 3, . . . , M — 1}. To determine the eigenvalues here one can proceed 
by employing well-known methods from polymer physics for the treatment of finite discrete 
chains lfT9l l28l. It turns out that H° also leads to eigenvalues of the form = 2 — 2 cos 9 X , 
but that now 9 X equals 9 X = irm/M, with m = 0, 1, . . . , M - 1 (25). 

These results allow us a straightforward analysis of the three network types given in Fig.[T]. 
Let us hence turn to considering the two-dimensional cases. 

3.1. The torus 

We start by considering a finite two-dimensional regular network with PBCs in both 
directions, i.e., a torus, see Fig.[TJc). In this case the Hamiltonian H* reads 

*F\jx,jy) = 2\jx,jy) - \ jx ~ l,jy) ~ \ jx + l,jy) 

+ 2\jx,jy) ~ \jx,jy - 1) - \jx,jy + 1), (14) 

where we take the labels j x and j y modulo M and N, respectively. Hence 

H t \j x ,j y )=H t (\j x )®\j y )) 

= (hsu.» ® \j y ) + \jx) ® (a p y \jy)) = m + h:;) ./,../,)• d5) 

as can be seen by directly inserting Eq. (fTTT) into Eq. (fT5t . by which Eq. (fl4l) is recovered. The 
result corresponds, of course, to a well-known result from solid-state-physics lfT71[T8l . Now, 
from H* = + H^ it follows that the eigenvalues of H* are given by 

Ae = \ dx + \ 9y = 4 - 2 cos 9 X - 2 cos 9 y (16) 
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with 9 X = 2mir/M with m = 0, 1, . . . , M - 1 and 9 y = 2nn/N with n = 0, 1, . . . , JV - 1, 
and that the eigenfunctions are 

M",JV 

\i, e ) = \M ® |^ v ) = 7^ E ex P H(0 (17) 



where ■ j = 9 x j x + By inverting Eq. (fl7l) . the state \ j) can be written as a linear 

combination of Bloch states \ipg). The projection of the eigenstates \ipg) on i.e., ipe(j) = 
(j\i>e) = exp[-i(0 ■ leads to ip e (j x + l,j y + 1) = expf-z'^a, + 9y)\^ {j xi jy), 

which corresponds to Bloch's theorem. The Bloch eigenstates of the system are orthonormal, 
i.e., (ipg\ip ') = 5 e e i, and they represent a complete set, i.e., J2e \'4'e)('4 , e\ = 1 holds. 
The quantum mechanical transition amplitude from state \j) to state \k) in time t now reads 

a k,j{t) = i JZ^e'l exp[-z0' ■ fc] exp(-iHt) exp[i0 • j] | Vo> 

0,0' 

= ^Y1 e M~^et) exp[-z6» • (k - j)}. (18) 
e 

For networks of infinite size, i.e., when M, N — > oo, the sums of Eq. (fT8t turn into integrals. 
This leads to 



7T 

M lmi^a fc j(t) = 6XP j^ 4 ^ y cf^ exp[-i0 a! (A; a . — exp(i2tcos9 x ) 

— 7T 

TT 

x / d0 y exp[— i9 y (ky — j y )\ exp(i2t cos9 y ) (19) 

— 7T 

= ik.-j.ik.-j, e ^{-iAt)J kx _ jx {2t)J ky _ jy {2t), 

where Ji(2t) is the Bessel function of the first kind ll29ll . In agreement with previous 
calculations ifTOl ITT1 IT2l the transition probability from j to k reads then 

lim 7r kd (t) = [J kx - Jx (2t)Jk y - jy (2t)} 2 - (20) 

MjJV — >oo 

Remarkably, the lower bound for the quantum mechanical probability to be still or again 
at the initial node [see Eq. (ITOb l becomes exact for PBCs QUI . Namely, inserting Eq. (fT8t for 
k = j into Eq. Q results in 

7f(*) = ^5ZS exp (~ iAet ) exp(iA /t) 
fe 0,0' 

= ^E ex PH(Afl - A*)t] = (21) 
0,0' 
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3.2. The cylinder and the rectangle 

We continue by considering now rectangular networks, on which we may or may not apply 
PBCs. In both cases for interior points Eq. (fl4b holds. Let us start from the cylinder, for which 
we impose PBCs in the x- and OBCs in the y-direction. Eq. (TBI) then also holds for all j x 
modulo M, where we identify \M,j y ) = \0,j y ) and \l,j y ) = \M + 1, j y ). For sites on the 
upper and lower row of the cylinder one has: 

H c \j x , 1) = 3\j x , 1) - \j x -1,1)- \j x + 1,1)- \j x , 2) (22) 

and 

H%, N) = 3\ Jx , N) - \j x — 1,N) — \ Jx + 1,N)- \j x , N — 1). (23) 
It is now a simple matter to show, in analogy to Eq. (fl~5l) that for all \j x ,j y ) one has 

ir. /,•./,} = (H£|j»» ® |j»> + ® (!«;;./„)) = (h* + 1 1;;). /,../,;)• (24) 

Hence, H c = + H° and the problem separates. Now the relation \ c e = + \ o 0y = 
4 — 2 cos 0-j - 2 cos ^ holds, where 6> x = 2nm/M with m = 0, 1, . . . , M - 1 and fty = 7rn/A^ 
withn = 0,1,..., JV- 1. 

Analogously, for a rectangle, we have OBCs both in the x- and in the y-direction. Let us 
denote the corresponding hamiltonian by H r . For H T the right-hand side of Eq. (TBI) holds for 
the internal nodes and the right-hand sides of Eqs. (l22l) and (l23l) hold (apart from the corners) 
for the upper and lower rows. Furthermore, similar expressions hold for the nodes on the left 
side and on the right side of the rectangle (again, excluding the corners). This means that for 
all nodes considered, we have 

W\j x ,j y ) = (H° x + K° y )\j x ,j y ). (25) 

Now it remains to be shown that Eq. d25t holds also for the corners. Exemplarily, we consider 
the corner, |1, 1), for which one has 

tT|l,l) = |M)-|2,l) + |l,l)-|l,2) 

= (H°|l)) ® |1) + |1) ® (H°|l)) = (H° + H°)|l, 1). (26) 

This completes the proof of Eq. d25t for all \j x ,j y ), which means that H r = H° + H° and 
thus the problem separates. Hence, the eigenvalues of H r are given by = A^ + A^ = 
4 — 2 cos 9 X — 2 cos 9 y , where 9 X = nm/M with m — 0, 1, . . . , M — 1 and 9 y = nn/N with 
n = 0,l,...,iV-l. 
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FigureEldisplays the eigenvalues for a rectangle, for a cylinder and for a torus containing each 
15 x 11 sites. Now, in each case the eigenvalues lie with in interval [0,8[, see e.g. Eq. (fT6l) . 
Figure |3 moreover, shows nicely that the n-th eigenvalue of the rectangle is always smaller 
than or equal to the n-th eigenvalue of the cylinder, which in turn is smaller than or equal to 
the n-th eigenvalue of the torus. 




Figure 2. Eigenvalues X n , arranged in ascending order for the 15x11 rectangle, cylinder and 
torus, respectively. 



4. Probability of being at the original site 

Classically, the average probability to be still or again at the initial site is given by p(t) 
[Eq. ©]. For our structures, we can get p(t) without numerically diagonalising A, because 
p(t) depends only on the eigenvalues, which are exactly known from Eq. (fT6l) (with different 
values for 9 X and 9 y , depending on the BCs of the structure). Quantum-mechanically the 
corresponding expression is 7f (t) [Eq. ©], to which [Eq. (fTOb l was shown to be a lower 
bound. We recall that 7f(t) depends also on the eigenvectors, whereas fx(t) does not. In the 
following figures we consider as an example the 15x11 network with the three different BCs. 
We start from the torus, i.e., from PBCs in both directions. Figure in which we plot p(t), 
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7f(£) and n(t), confirms the fact that for PBCs n(t) and /j,(t) coincide [see Eq. (ETTl l. In the 
intermediate range (in Fig.|3]from t = 0.5 to t = 5) the classical probability p(t) is algebraic, 
i.e., here we have p(t) ~ t" 1 . 




t 



Figure 3. Classical, pk,k(t)> an d quantum mechanical, TTk,k(t), probabilities to be still or 
again at the initial site and the lower bound ji(t) defined by Eq. (II Oi for the 15 x 11 torus. 




Figure 4. Classical, pk,k(t), and quantum mechanical, 7Tfe,fe(£), probabilities to be still or 
again at the initial site and the lower bound fi(t) defined by Eq. d!0> for the 15 x 11 cylinder. 

Figures |4] and |5] show for the cylinder and for the rectangle, respectively, the probability to 
be still or again at the initial site in the classical and in the quantum case. Furthermore, we 
also show the quantum mechanical lower bound /i(t), which now differs from 7f(t) . However, 
depending on the BCs, this difference is very small and the maxima and the minima of fj,(t) 
and 7f(£) occur basically at the same times. By comparing Fig.|3]to Fig.HJand to Fig.|5]one 
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Figure 5. Classical, pk,k(t), and quantum mechanical, 7Tfe,fe(£), probabilities to be still or 
again at the initial site and the lower bound /i(t) defined by Eq. J10I for the 15 x 11 rectangle. 



can conclude that pit) reaches the asymptotical behaviour given by pit) = 1/Af earlier for 
the torus than for the cylinder and for the rectangle. We explain this finding by the fact that 
some sites on the torus are more connected than the corresponding sites on the cylinder or on 
the rectangle. 

5. Limiting probability distributions 

The unitary time evolution prevents the quantum mechanical transition probability from 
having a definite limit when t — > oo. For comparison to the classical long time probability, 
the long time average of n k j, i.e., 

1 <- T 



Xkj = Jim - / dt TT kJ (t), (27) 

1 — >oo ± 







can be used, see for instance [31 1. In order to avoid the numerical integration, which requires 
long computing times already for moderately large network sizes, we rewrite the LP as follows 

G3 



T 



1 f 

Yfc ,• = lim — dt 

A ,3 T^ooT J 



^2(k\exp(-iH.t)\q n )(q n \j) 

n 

= J2^ q n)( ( ln\j)(j\ ( l m )( ( lrn\ k ) I ^ \ J dt exp[-l{\ n - \ m )t] 

= ^ S ^(k\q n )(q n \j)(j\qJ(q m \k). (28) 

n,m 

This expression simplifies the numerical evolution considerably. 
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Let us first consider the behaviour of Xk,j taking as initial node j the centre of a M x M 
network, where M is odd. The classical behaviour is diffusive and it leads to an equipartition 
among all nodes of the network, hence classically lim t _^ 00 p(t) = 1/jV. Quantum- 
mechanically however, the situation is different. Figure |6] shows the pattern of Xk,j f° r a 
15 x 15 cylinder. There is a very pronounced global maximum at the centre. Other local 
maxima are to be found on the nodes placed on the diagonals of the structure and on the 
shortest paths which connect the centre to the four sides. The LPs of the other nodes produce 
a plateau which is approximately half of the local maximum. Furthermore, we find the same 
result, up to our numerical precision, when we consider a square or a torus with the same 
initial condition. 



Xk,centre 



0.035 
0.025 
0.015 
0.005 




Figure 6. Limiting probabilities for a 15 x 15 cylinder, the walk starts at a central node. 



Now we place the initial excitation at one corner of the structure, say at j = (1,1). For the 
torus all nodes are topologically equivalent, so that the result is identical to the one displayed 
in Fig.0 We expect hence differences only in the case of squares and of cylinders. 



Quantum transport on two-dimensional regular graphs 14 

Now, for square networks one intuitively expects the Xk,j to be symmetric with respect to 
the centre of the network, i.e., one expects the node (j x ,j y ) and its mirror node, defined by 
(M + 1 — j x , M + 1 — j y ) to have the same LP. In general, this turns out to be true, but, 
there are some exceptions lfT2ll . namely for M — 6, 12, 15, 18, 21, 24, 30, 36, .... In Fig.EJwe 
exemplify this behaviour. Thus the 14 x 14 square has a symmetric LP, shown in Fig. Eta), 
whereas the 15 x 15 square has an asymmetric LP, see Fig.Etb); note the differences along 
the borders and also the differences at the points (1, 1) and (M, M) in Fig.Etb), compared to 
Fig. Eta). These differences will become clearer in the discussion of Fig. [8] 






Figure 7. LP distributions for walks starting at j = (1, 1) for (a) the 14 x 14 square, (b) the 
15 x 15 square, (c) the 14 x 14 cylinder and (d) the 15 x 15 cylinder. 



With cylindrical boundary conditions we find a similar situation. However, the symmetry is 
of another kind. We choose OBCs in the y-direction and PBCs in the x-direction. In Fig. Etc) 
we plot the LPs for the 14 x 14 cylinder. In this case there is a mirror symmetry on the rows 
with constant k x , i.e., Xk,c is the same for the nodes (k x , k y ) and (k x , M — k y ). This is not 
true for the 15 x 15 cylinder, for which the LPs are given in Fig.Etd). 

Because of the PBC along the ^-direction one expects the same LPs, for instance, at the 
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second and the last nodes along the x-direction. This is indeed borne out by the figures. 
To render clearer these asymmetric behaviours, we plot in Fig. Ufa) as a function of M 
the differences Xc,c — Xoc,c f° r c = (1, 1) an d oc = (M,M). We obtain values different 
from for M = 6, 12, 15, 18, 21, 24, 30, . . .. The situation for the cylinders is shown in 
Fig. [Mb), where we plot Xc,c — X^y,c with c y = (1, M). The asymmetries appear then for 
M = 6, 15, 18, 21, 30, . . .. Remarkably, in both cases, non-zero Xc,c ~ Xcy,c values are found 
only for iV values which are multiples of three. 



0.04 



o 



0.03 - 



0.02 



0.01 




4 6 



10 12 14 16 18 20 22 24 26 28 30 
M 



0.04 



0.03 



0.02 



0.01 - 




4 6 



10 12 14 16 18 20 22 24 26 28 30 
M 



Figure 8. (a) Differences between the LP for CTQWs that start at c = (1, 1) to be at c, Xc,c, 
and to be at oc = (M, M), Xoc,c for OBC. (b) Differences between the LPs for CTQWs on a 
cylinder that start at c = (1, 1) to be at c, and to be at c y = (1, M). 



One can now sum the LPs along one of the two directions (x or y) of the network. For a 
cylinder this leads after summing along the open y-direction to a ring-like behaviour and after 
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12 14 16 



0.14 



1 0.12 




14 16 



Figure 9. Comparison between the marginal LPs to J^k Xk,c (full line) and J^k Xfc,c 
(dashed line) on a 14 x 14 network (a) and on a 15 x 15 network (b), see text for details. 



summing along the x-direction with PBCs to a chain-like behaviour. Figure |9t a) shows the 
LP sums Y^k x =i Xk,c (full line) and Ylk =i Xk,c (dashed line) for the 14 x 14 network with 
k = (k x , k y ) and c = (1, 1). Both curves have two maxima. For the chain the maxima occur 
at the extremes (k y = 1 and k y = 14), whereas for the ring they arise at k x = 1 and k x = 8. 
In Fig. |9tb) we plot the corresponding results for the 15 x 15 network. For the chain one 
finds two equal maxima at the extremal points and a low, equal probability in the middle. 
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For the ring, on the other hand, we find just one maximum at k x = 1 whose value is 
(2M - 1/M 2 ) « 0.1288 and a plateau (M - 1)/M 2 « 0.06222 for fc a ^ 1. Here, distinct 
from Fig.|9ta), the plateau has the same value in both cases (chain or ring). 

5.2. M x N networks 

We extend now the analysis of section 15.11 to general M x N networks with M ^ N, 
considering again different boundary conditions. Figure [TO] shows examples of symmetric 
and asymmetric patterns for networks of sizes 15 x 10 [(a) and (b)], 15 x 11 [(c) and (d)] 
and 16 x 10 [(e) and (f)], the networks being either rectangles [(a), (c) and (e)] or cylinders 
[(b), (d) and (f)]. For rectangles, OBCs, the global maximum is at the initial site and there are 
pronounced maxima of the LPs along the edges of the structure; this is similar to the situation 
for M x M-networks shown in Fig. IT] However, here we do not find local maxima along the 
diagonals. Clearly, the LPs for the 15 x 10 rectangle show more structure fFig. flOt a)! than 
the other two networks. Thus, the 15 x 11 network fFig. fToT c^l is symmetric and its LPs take 
only three different values; the 16 x 10 network fFig. flOt e)! is very regular and its LPs are 
concentrated on the peripheral nodes. 

Figures fTOt b). (d) and (f) present the corresponding cylindrical networks. The LPs for the 
15 x 10 network have their global maximum at c = (1,1) and, again, there is no mirror 
symmetry for the nodes (k x , k y ) and (k x , N—k y ). The Xk,c for the 15 x 11 cylinder rFig. flOl' d')! 
are very regular. One observes two equally large peaks at c = (1,1) and at c y = (1, 11). For 
the other points of the row k x — 1, the LPs are equal. The same holds for the nodes of the 
cylinder (at k y = 1 and at k y = 11) and for the other remaining nodes of the network. Figure 
IToT f) displays Xk,c f° r the 16 x 10 cylinder; the LP distribution is quite different from the 
previous two patterns. As in Fig. He) there appear two "crests", at k x = 1 and at k x = 8, a 
result due to the PBCs along the x-axis and to the fact that here M is even (M = 16). 
In order to systematically analyse the asymmetries of Xk,j in rectangular networks, we now 
fix N to be N = 15 and vary M, taking 4 < M < 30. Figure [TTT a) shows for rectangles the 
difference between the LPs on the initial corner Xc,c and on the corner Xoc,c- In the analysed 
range, 4 < M < 30, the value Xc,c — Xoc,c displays varying patterns. For 4 < M < 13 one 
can associate Xc,c — Xoc,c = to odd values of M and Xc,c — Xoc,c ^ to even values of 
M; for larger M the situation becomes more complex. Figure fTTT' b') displays the situation for 
cylinders with N = 15 and 4 < M < 30. Here we plot Xc,c — Xcy,c, with c y = (1, N). Figure 
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rnT c) shows the situation for cylinders with M =15 and 4 < iV < 30. This last case looks 
quite regular, with non-vanishing values only for iV = 10, 15, and 30. Hence, Figs. \T\\ b) 
and nTT c) (N < M) show for cylinders that changes of the radius lead to more asymmetric 
situations than changes of the length. 
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Figure 10. LPs for (a) the 15 x 10 rectangle, (b) the 15 x 10 cylinder, (c) the 15x11 rectangle, 
(d) the 15 x 11 cylinder,(e) the 16 x 10 rectangle, (f) the 16 x 10 cylinder. 
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Figure 11. (a) Rectangles M X N with N — 15 and varying M: Differences between the 
LPs for CTQWs that start at c = (1, 1) to be at c, Xc,c, and to be at oc = (M, N), Xoc,c- (b) 
Cylinders M x N with N = 15 and varying M: Differences between the LPs for CTQWs 
that start at c = (1, 1) to be at c, and to be at the node c y = (1,N). (c) Cylinders M x N 
with M = 15 and varying N: Differences between the LPs for CTQWs that start at c = (1,1) 
to be at c, and to be at the node c y = (1, N). 
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We close by briefly reviewing some of the details of the calculations of the LPs, which give an 
indication of the origin of the asymmetries. As we have seen in Sec. the Hamiltonian of the 
problem separates in the two directions. Using this fact, it is easy to show that the transition 
probabilities 7Tk,j{t) can be written as the product of the two separate probabilities for each 
direction, i.e. 



with a kx ,j x {t) = Y^e x e *v{-i^eJ){K\$o x ){^o x \jx) , and similarly for the y-direction. 
Now, according to Eq. (l28l) . the LPs are given by 

i r T 

Xkj = Jim - / dt ir kxJx (t) ir ku j y (t) 



where F^j is a time independent function, which depends on the eigenstates associated 
with 9 X , 6' x , 9y, and 0' y . Because of the limit in the time integral in Eq. d3*0"l) there are only 
contributions to Xkj if a value (X 6x — \$> x ) for the s-direction has a counterpart — (X 6y — A#/ ) 
in the y-direction. 

A careful analysis of the differences (Xq x — Xg x ) indicates where the asymmetries stem from. 
For finite chains we obtain (Xq x — Xq/) = 2 cos 6^ — 2 cos 8 ] x . For simplicity we consider 
now finite M x M networks with OBCs, see |[T2|, because then the eigenvalues are the 
same in both directions. It turns out that for 6 X ^ 6' x the value (A^ — Xg x ) appears only 
once or twice for all symmetric cases. However, for the asymmetric cases, some of the 
(A^ — Xg x ) values (again for 9 X ^ 9' x ) appear more than twice. Therefore, there are more 
contributions to Xkj in the asymmetric cases than in the symmetric cases. Given that we 
have all eigenvalues analytically, we can in principle study the asymmetries fully analytically. 
However, a complete analysis, including different boundary conditions, requires a thorough 
investigation of all possible differences of eigenvalues. This is a quite extensive task and is 
clearly beyond the scope of this paper. Thus, a detailed analysis of the asymmetries will be 
given elsewhere. 

6. Conclusion 



(29) 




To summarise, we have systematically analysed the role of different boundary conditions on 
the quantum mechanical transport on two dimensional graphs. We have kept the systems 
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as simple as necessary to still highlight the complex behaviour of the quantum mechanical 
transport. In the long time average of the transition probability distribution Xk,j we found 
asymmetries, depending on the size and the BCs of the graph. Although graphs with toroidal 
topologies have a perfectly symmetric Xk,j, this does not have to be the case for graphs of the 
same size with open BCs. For networks with an equal number of nodes in both directions 
as well as for networks with an unequal number, the asymmetries are best visualised by 
comparing the LP of the initial node to the one of an appropriate "mirror node". Here we 
found significant differences for the specific sizes of the network depending on the BCs. We 
further showed that summing along one side of the graph leads to the known behaviours for 
quantum mechanical transport over chains (OBCs) or rings (PBCs). 
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